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Abstract 

The central smoothness of the functions defining a LTB solution plays a crucial role in their ability to 
mimick the effects of the cosmological constant. So far attention has been focused on models while 
in this paper we approach it a more general way, investigating the implications of higher order central 
smoothness conditions for LTB models reproducing the luminosity distance of a AC DM Universe. Our 
analysis is based on a low red-shift expansion, and extends previous investigations by including also 
the constraint coming from the age of the Universe and re-expressing the equations for the solution 
of the inversion problem in a manifestly dimensionless form which makes evident the freedom to 
accommodate any value of Hq as well, correcting some wrong claims that the observed value of Hq 
would be enough to rule out LTB models. 

Higher order smoothness conditions strongly limit the number of possible solutions respect to the 
first order condition. Neither a or a C* LTB model can both satisfy the age constraint and mimick 
the cosmological constant for the luminosity distance. One difference is in the case in which the age 
constraint is not included and the bang function is zero, in which there is a unique solution for 
models but no solution for the C* case. Another difference is in the case in which the age constraint 
is not included and the bang function is not zero, in which the solution is undetermined for both 
and C* models, but the latter ones have much less residual parametric freedom. 

Our results imply that any LTB model able to fit luminosity distance data and satisfy the age con- 
straint is either not mimicking exactly the AC DM red-shift space observables theoretical predictions 
or it is not C°° smooth. 



*Electronic address: |aer@phys.ntu.edu.tw| 
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I. INTRODUCTION 

A vast range of cosmological observables jl-^ when interpreted under the assumption of 
large scale spatial homogeneity supports a dominant dark energy component, giving rise to a 
positive cosmological acceleration. This value of the cosmological constant is different from 
the one predicted by quantum field theory by few orders of magnitude. Given this evident 
inconsistency about our understanding of the Universe, a lot of attention has been devoted to 
alternatives to dark energy, such as modified gravity theories or large scale inhomogeneities. 
These cosmological models are based on violating some of the assumptions of standard cosmol- 
ogy. One of these assumptions is that the Universe is well described by general relativity at 
any scale, and a violation of this leads to modified gravity theories, which consist of theories of 
gravity which agree with general relativity on small scales such as the solar system, and depart 
from it on larger cosmological scales. The other main assumption of the standard cosmological 
model is that the Universe is homogeneous on sufficiently large scales, and can be correctly 
described by a spatially homogeneous and isotropic metric, i.e. the FLRW metric. Since there 
is no conclusive evidence of such a large scale homogeneity, and the Copernican principle is 
more of a philosophical principle than an experimental evidence, inhomogeneous cosmological 
models have been studied as alternatives to dark energy. 

It has been proposed j^, [lO| that we may be at the center of an inhomogeneous isotropic 
universe without cosmological constant described by a Lemaitre-Tolman-Bondi (LTB) solution 
of Einstein's field equations, where spatial averaging over one expanding and one contracting 
region is producing a positive averaged acceleration a^, but it has been shown how spatial 
averaging can give rise to averaged quantities which are not observable ll| . A method to map 
luminosity distance as a function of redshift Dl{z) to LTB models has been proposed 13|, Il5| . 
showing that an inversion method can be applied successfully to reproduce the observed Dl{z). 
Analysis of observational data in inhomogeneous models without dark energy and of other 
theoretically related problems is given for example in 

In this paper we study the anaytical solution of the inversion problem (IP) which consists in 
matching exactly the terms of the low redshift expansion for the relevant cosmological observ- 
ables, in particular the role played by central smoothness. As observed in 27|] the smoothness of 
the inhomogeneity profile is important to allow to distinguish these LTB models from ACDM, 
which is consistent with the conclusion that the solution of the IP is possible for both the 
luminosity distance Dl{z) and the red-shift spherical mass energy p(^z) only if we allow a not 
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smooth radial matter profile [28| , while smooth models 29| can be distinguished from ACDM, 
at least in principle, without any ambiguity. In these previous investigations the smoothness 
condition consisted in assuming that the function necessary to define the LTB geometry did 
not contain any first order derivative term at the center. In this paper we show that in order to 
ensure smoothness at higher orders it is necessary to impose the condition that any odd deriva- 
tive of the functions defining the model vanishes at the center, and study its implication on 
the solution of the inversion problem for Dl{z), including for the first time, also the constraint 
coming from the age of the Universe. 

It is important to observe that according to our definition of the IP we look for the LTB 
models which match the coefficients of the redshift expansion of the observables corresponding 
to the best fit flat ACDM models, but since these models depend only on the two independent 
parameters Hq, ^m, LTB models have a higher number of parameters, implying that they could 
actually provide an even better fit of experimental data. Our definition of IP is nevertheless 
quite natural from a mathematical point of view since it consists in matching the red-shift space 
theoretical predictions of different models within the range of validity of the Taylor expansion/. 
The analytical formulae can also be used for low-redshift data fitting, with the advantage of not 
depending on any functional ansatz. We can for example draw contour plots for the coefficients 
of the expansion of the functions defining the LTB model , in the same way we do for As far 
as observations are concerned it would in be more important to fit directly actual experimental 
data rather than mimic the best fit theoretical ACDM model, so our conclusions should be 
considered keeping this in mind, and we leave to a future work the analysis of experimental 
data. 



This kind of analysis has been performed by different groups 3(J ? , |3l|, and recently 32 1 
it has been claimed that accurate measurements of Hq are sufficient to rule out best fit void 
models. This claim has been corrected in jsSi, and the dimensionless form in which we re- write 
the equations to solve the IP, makes it evident even in the low-redshift expansion approach we 
adopt, that Hq^^ can accommodate any observed Hq^^. 

The paper is organized as follows. We devote a section to the general requirements for the 
LTB solution to be smooth, i.e. infinitely differentiable, at the center, generalizing it to the 
requirement that all the odd derivatives of the functions of the radial coordinates defining the 
model should vanish at the center; without this smoothness condition even derivatives of the 
energy density would diverge at the center. Using a local Taylor expansion we discuss in details 
the existence, uniqueness and degeneracy of the solutions of the IP consisting in matching the 



observations of the luminosity distance, Hq and the age of the Universe. We classify models 
according to the central smoothness and the presence of an isotropic but inhomogeneous big 
bang. We perform the some classification using only the first order smoothness condition , and 
the the higher order ones, showing how these latter strongly impact the number of smooth 
solutions. Higher order smoothness conditions in fact strongly limit the number of possible 
solutions respect to the first order condition. One difference is in the case in which the age 
constrain is not included and the bang function is zero, in which there is a unique solution 
for models but no solution for the C* case. Another difference is in the case in which 
the age constraint is not included and the bang function is not zero, in which the solution is 
undetermined for both and C* models, but the latter ones have much less residual parametric 
freedom. This implies that any other LTB model able to fit luminosity distance data, satisfy 
the age constraint and fit some other observable is either not mimicking exactly the ACDM 
red-shift space observables theoretical predictions or it is not C°° smooth. 



II. LEMAITRE-TOLMAN-BONDI (LTB) SOLUTION 



This is one form to write the Lemaitre-Tolman-Bondi 
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where i? is a function of the time coordinate t and the radial coordinate r, R = R{t,r), E is 
an arbitrary function of r and R,r = OR/ dr. 
from the Einstein's equations we obtain 
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with M = M{r) being another arbitrary function of r and the dot denoting the partial derivative 
with respect to t, R = dR{t,r)/dt. It is possible to introduce the variables 
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SO that Eq. ([T]) and the Einstein equations (E]) and (El) are written in a form similar to those 
for FLRW models, 
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The solution takes now the form 
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where i] = rr = /* dt' /A(t', r) and A(t{r], r), r) = a{r], r) and tb(r) is another arbitrary function 
of r, called the bang function, which is an integration constant corresponding to the fact that the 
big-bang initial singularity can happen at different times at different distances from the center. 
The inhomogeneity of the location of the singularities can lead to possible causal separation 
ll| between the central observer and the spatially averaged region for models with positive 
an- There is an open debate about the presence of the bang function in LTB models consistent 
with observations, since the bang function would correspond jl2, 13] to large decaying modes in 



the primordial curvature perturbations in disagreement with the cosmic microwave background 
(CMB) observations. Other authors jl^ claim instead that the present CMB observations are 
not incompatible with an isotropic but inhomogeneous bang function. In this paper we will 
not deal with the physical implications of an inhomogeneous big bang, but in the section in 
which we analyze the possible solutions to the IP we will distinguish between the homogeneous 
and inhomogeneous big bang case in order to provide a general analysis of all the possible 
mathematical cases. In the rest of paper we will use this last set of equations. Without loss 
of generality we can set the function po{r) to be a constant, po{r) = po = constant, which the 
choice of coordinates in which M(r) oc r^. 



III. THE INVERSION PROBLEM FOR 



The general low re-shift formula for the luminosity distance for a central observer is 28] 
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It is convenient to re-write it in terms of dimensionlcss parameters according to: 
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to get for example: 
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As expected observable quantities such as Dl{z) do not depend on oq. Based on the above 
results it is convenient to schematically write the coefficients of the expansion of D^^[z) in 
the form: 
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which clearly shows what is the dependence between the coefficients of the expansion of the 
function defining the LTB model K{r),Th{r) and those of the red-shift space observable. 

The analytical version of the inversion problem up to order n consists in matching the 
coefficients of the expansion of D^^^^'^[z) and D^^[z): 

^LTB ^ ^KCDM ^ l<^<n. (22) 

IV. CENTRAL SMOOTHNESS AND OBSERVER LOCATION FINE-TUNING 



It has been proved using a local Taylor expansion approach that 29| only not centrally 
smooth LTB models can mimick the effects of the cosmological constant for more than one 
observable. This is independent from the value of the cosmological constant, and is related 
to the number of free independent parameters on which the functions defining the solutions 
depend on. Central smoothness of the LTB solution is also important because is related to the 



sign of the apparent value 42| of the deceleration parameter q'o'^ . For this reason we will later 
classify the solutions of the inversion problem according to central smoothness, since the space 
of not smooth solutions is larger than that of smooth solutions, and in same cases there exist 
no smooth solutions, as shown in the next section. 

From a general point of view a function of the radial coordinate /(r) is smooth at the center 
r = only if all its odd derivatives vanish there. This can be shown easily by looking at the 
the partial derivatives of even order of this type for example: 



afafaf/(v/a;2 + 2/2 + ^2)^ (23) 

where {x,y, z} are the cartesian coordinates related to r by = + + z"^. Quantities of 
the type above diverge at the center if d^'^~^^f{r) ^ for 2m + 1 < 2n. If for example the ffist 
derivative f'{0) is not zero, then the laplacian will diverge. 

The general central smoothness conditions for a function of class C* are: 



k2m+l = 0, (24) 
^2"^+! = 0, (25) 
2m + 1 < z, (26) 

which must be satisfied for all the relevant odd powers coefficients of the central Taylor expan- 
sion. For definiteness we will present here the results of the calculations to the second order 
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in r, corresponding to third order in z for Dl{z). We will show later nevertheless that our 
conclusions about the existence and uniqueness of the solution are independent of the order at 
which we truncate the above expansion, and that in general the expansion of the Dl{z) to the 
i-th order in z requires the expansion to (i-l)-th in r for k{r) and tfe(r). 

A possible criticism about the physical relevance of considering central smoothness as a cri- 
terion to study LTB models is that the available observations, such as Supernovae for example, 
are at redshift higher than 0.013, i.e. relatively far from the central observer. This implies 
that fitting data with centrally smooth models can still give good results since there are not 
sufficiently low redshift observations to affect statistically the data fitting. If we assume nev- 
ertheless that the ACDM model is also phenomenologically valid at redshift lower than the 
range of the currently available cosmological observations, we should expect that the goodness 
of the fit at very low redshift should be quite poor for smooth LTB models, as predicted by the 
low redshift expansion analysis. This is the reason why smooth models can give good results in 
fitting supernovae data, even if their low red-shift luminosity distance is not mimicking exactly 
the ACDM model, since most of the data used to compute the goodness of the fitting are at 
higher red-shift where the central smoothness is not so important anymore. On the other side 
the existence of a center is a distinctive feature of LTB models and consequently the behavior 
of the solution around the center cannot be neglected. 

From a mathematical point of view it has been proved 29| that central smoothness is a key 
characteristic for a LTB model to be able to mimick ACDM for more than one observable, 
independently from the value of the cosmological constant. From a data fitting point of view, 
this cannot be checked yet due to the lack of sufficiently low redshift observations, which does 
not diminish the relevance of such an important mathematical feature of LTB models as viable 
cosmological models. 

Finally we note that the LTB models which have been claimed to be incompatible with 
recent accurate Hq measurements are based on the assumption of the central location of the 
observer, which we also make in our analysis. In these models the observed dipole anisotropy 
of the CMB radiation is associated to our velocity relative to the CMB, and is the same as 
in a ACDM model, since for a central observer there is no additional anisotropy added by 
the propagation of photons in a spherically symmetric inhomogeneous space. Alternatively 
the dipole anisotropy component of the CMB could be considered the effect of an off-center 
position of the observer in a spherically symmetric space, due to the fact that space is not 
isotropic around the off-center observer. Combining Supernovae observations and the CMB 
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anisotropy dipole it has been found 16| that the location of the observer has to be within one 
percent of the void radius, a fine tuning which is in clear violation of the Copernican principle, 
and remains the main criticism against these models. Independently from the position of the 
observer LTB solutions have a center, and the smoothness conditions would still apply even to 
the case of an off-center observer. While it could be interesting to investigate how an off-center 
position of the observer could allow to improve the solution of the inversion problem, it is clear 
that since the only additional parameter would be the distance from the center, the general 
conclusions achieved in this paper will not be substantially affected up to an additional residual 
degree of freedom in determining the solution. 

V. HOW MANY INDEPENDENT PARAMETERS DETERMINE LOCALLY A LTB 
DUST MODEL? 

In order to fit observation with a LTB dust model we need to understand well how many 
are the really independent free parameters. From the previous section we know that that any 
expression containing po, rjo, can be re-expressed in terms of ag, Hq, go, and that introducing 
appropriate dimensionless quantities does not appear explicitly in observables such as Dl{z). 
This suggests that there are only two independent parameters to fix in order to fit observations, 
and we will prove this in a more rigorous way. 

We can start from rewriting the equations eq. (11111131) in the following form: 



1 = ^^^ = fi{qo,KoM, (27) 

Ko go 



1 

^~ 2^ 



1 = ^ (^^) ^ = /2(go, Ko, Qm) , (28) 

2 



-- h{qo,Ko,n), (29) 
where we have used introduced the dimensionless parameter Qrn and Tq according to: 



Po = SalH^Qm , (30) 
r]o = ToiaoHo)-'. (31) 

Another possible constraint could come from the age of the Universe. Since the dust models we 
are considering cannot be used to describe the universe during the radiation dominated stage, 
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they cannot really be considered fully realistic cosmological models. Nevertheless, since most of 
time the Universe has been matter dominated, wc can safely neglect the contribution to the age 
of the Universe coming from the radiation dominated era. Using the exact analytical solution 
t{r], r) we can write: 



Tu^Hot{rio,0)^^ 

ZKq 



1 . 



^U{n^,Ko,To), (32) 



where we have introduced the dimensionless quantity Tu = tuHo, which is the age of the 
Universe in units of Hq^. 

We have four equations four parameters Ko,To,qo, Jl^, while Hq remains free. This is man- 
ifest from the fact that the equations are written in terms of dimensionless parameters which 
are independent from the value of Hq. 



VI. Ho OBSERVATIONS AND LTB MODELS 

The value of Hq is defined observationally in a model independent way as: 

(Ht)-' = lim . (33) 

^ " ^ z^O dx x=z 

Let us now expand locally the observed luminosity distance : 



Dt^iz) = Df^z + Dtz' + .. = j^^z + Dt + .. (34) 

In order to avoid confusion we will denote here all the quantities related to a LTB model with 
superscripts according to: 



H^^''{n,^,Ko,To), (35) 

(36) 

We will also limit ourselves to the homogeneous bang function case, i.e. set ti = 0, and mention 
later how including the bang function will make our arguments even stronger. According to 
the expansion for the luminosity distance obtained previously we have 
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1 ~ jj-LTB ' (2''') 

HLtBj^ltb ^ Uq^^^,K,X), (38) 
HLTBj^ltb ^ f,(g^^^,K,,K„T^,T^), (39) 

(40) 

The advantage to write in this form the coefficients of the expansion for the luminosity distance 
is that in this way we can clearly understand the role of Hq as the natural scale of the problem, 
which is set by Hq^^. 

Contrary to some previous attempts l?! to fit observation with void models the correct 
way to interpret this equations is not to look for the particular set of parameters which give 
Hq^^ = Hq^^, but to directly fix Hq^^ = H^^'^ and then look for the dimensionless parameters 
which solve the inversion problem, i.e. specify the LTB model in agreement with observations. 
This means that there is always a set of LTB solutions in agreement with any value of H^'^'^ , as 
long as low redshift observations for the luminosity distance are concerned, and an arbitrarily 



precise measurement of Hq^'^ is not by itself able to rule out void models[32[], but only able to 
establish accurately the natural scale in terms of which to define the model. 

In order to design a model in good agreement at low redshift with observations we can match 
the Taylor expansion of the observed luminosity distance with the corresponding LTB terms : 



Combining the above equations with the ones derived in the previous section we finally get: 
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(49) 



As mentioned previously the advantage of writing the equations in this form is to clearly identify 
the role of Hq^^ as the scale in terms of which to define the LTB model, which is given directly 
from the first equation. We will analyze in the rest of the paper the different possible cases 
according to central smoothness and the presence of a bang function tf,(r). One last important 
point is that Hq^^ is a purely local quantity, i.e. only depends on the central value of k{r),ri 
and Qm, implying that the Taylor expansion approach we adopted gives exact results as long 
as Hq is concerned. 

VII. INVERSION PROBLEM SOLUTIONS FOR MODELS 

The number of equations to be satisfied to ensure the matching of the luminosity distance up 
to i-th order is 4 + (i — 1), where we have (« — 1) constraints because we consider the matching 
of the first term Df"^^ equivalent to fixing the scale in terms of which to define the rest of the 
parameters of the model. 



Age constraint 


LTB type 


Not C^, h^O 


Not C^, th = 




C\ h = 


Free parameters 


4+2(i-l) 


4+(i-l) 


4+2(i-l)-2 


4+(i-l)-l 


Constraints 


3+(i-l)+l 


3+(i-l)+l 


3+(i-l)+l 


3+(i-l)+l 


Solution 


Undetermined 


Unique 


No 


No 



TABLE I: The number of parameters on which LTB models depend expanding the luminosity distance 
to i-th order are shown according to the central behavior of the energy density (smooth/not smooth) 
and the presence of a homogeneous bang function tfc(r). In this case we are including the constraint 
coming from the age of the Universe. 

In order to understand the content of the tables it can be seen from eq. fl43ti45|) that in- 
dependently of the order of the matching of the Taylor expansion there always are four free 
parameters qo,KQ,TQ, present and three constraints. When including the age of the universe 
then we have an extra constraint. From the matching of the coefficients up to i-th order of the 
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No Age constraint 


LTB type 


Not C\ tb^O 


Not tb = 


C\ tb / 


C\ tb = 


Free parameters 


4+2(i-l) 


4+(i-l) 


4+2(i-l)-2 


4+(i-l)-l 


Constraints 


3+(i-l) 


3+(i-l) 


3+(i-l) 


3+(i-l) 


Solution 


Undetermined 


Undetermined 


Undetermined 


Unique 



TABLE II: The number of parameters on which LTB models depend expanding the luminosity distance 
to i-th order are shown according to the central behavior of the energy density (smooth/not smooth) 
and the presence of a homogeneous bang function t(,(r). In this case we are not including the constrain 
coming from the age of the Universe. 

luminosity distance expansion instead we get 2{i — 1) parameters corresponding to Ki_i,Ti_i, 
coming from Df^^ . Imposing the condition or setting to zero the bang function tb(r) reduces 
the 2{i — 1) parameters to the values reported in the tables. 

The condition of the existence of a solution of the inversion problem is: 

j^parameters > ^constraints , i > 2 (50) 

In the case of equality the solution is unique. 

A. Analysis of the relevant cases 

In the following it should be noted, as mentioned above, that we are assuming by construction 
that Hq^^ = (which is always possible), and that all the other parameters defining a LTB 
are dimensionless. The existence of a solution of the inversion problem can be easily determined 
by looking at the number of free parameters and the number of unknowns reported in the tables 
but there are some cases in which some extra explanations are required. 

1. with age constraint 

When the functions defining the model are at the center, i.e. when Ki = = 0, we have 
no solution if we impose the constraint coming from the age of the Universe, with or without 
the bang function tb{r). Apparently for i > 3 we have more free parameters than constraints, 
but for i = 2 we have that the system of five equations (143M7P has only four unknowns and 
the higher order parameters Ki_i,T^_^ appearing in Df^^ do not have any effect on equations 
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( I43ll47l) . If we do not include the constrain from the age of the Universe in eq.( l46l) . the solution 
is actually undetermined, because in this case the system of four equations ( l43ll45K71) has a 
unique solutions, but higher order terms will introduce extra degeneracy, since every term Di 
will contain two independent parameters Ki_i,Tf'_i, and only one more constraint equation. 

2. Not without bang function and with age constraint 

The solution is unique in this case, since the number of parameters is equal to number of 
constraints at any order. This implies that any other LTB which which can fit observational 
data is either not mimicking exactly the the cosmological constant or, if it does, it is not smooth. 
This is due to fact that we are imposing very restrictive constraints on the functions defining 
LTB models, while data fitting allows for much more freedom due to the tolerance introduced 
in the statistical analysis, which can allow good data fitting even if the low-redshift observable 
are not exactly the same, and also because there are very few data at very low-redshift, making 
this mismatch not statistically important in evaluating the goodness of the fit. 

3. Not case with bang function 

When we allow for not models with bang function, we gain some local extra parametric 
freedom , which implies that the age of the Universe and luminosity distance are not enough 
to uniquely determine a LTB model in agreement with observations, leaving more freedom to 
fit other observables. 

VIII. INVERSION PROBLEM SOLUTIONS FOR C' MODELS 

In this section we will consider the implication of the C* conditions on the solution of 
the inversion problem. In general the space of free parameters is reduced further respect to 
the case, making the solution of the inversion problem some time impossible. As already 
observed for the models it can be seen from eq. fl43H50l) that there are four free parameters 
qo,Ko,To, Qrn present independently of the order of the matching of the Taylor expansion, 
and 2{i — 1) parameters corresponding to i^j_i,Tj_i, coming from Dj^^^. Imposing central 
smoothness or setting to zero the bang function tf,(r) reduces the 2{i — 1) parameters to the 
values reported in the tables. The number of constraints are the same as the ones in the 
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previous sections. 



Not smooth 


Age constraint 


Yes 


No 


LTB type 




tb = 




tb = 


Free parameters 


4+2(i-l) 


4+(i-l) 


4+2(i-l) 


4+(i-l) 


Constraints 


3+(i-l)+l 


3+(i-l)+l 


3+(i-l) 


3+(i-l) 


Solution 


Undetermined 


Unique 


Undetermined 


Undetermined 



TABLE III: The number of parameters on which LTB models depend expanding the luminosity dis- 
tance to i-th order are shown in the not centrally smooth case according to the presence of a homo- 
geneous bang function tb{r) and the age constraint. 



Smooth without age constraint 


Matching to i-th order in Dl{z) 


Even 


Odd 


LTB type 


Smooth, tb ^0 


Smooth, tb = 


Smooth, tb ^0 


Smooth, tb = 


Free parameters 


4+(i-2) 


4+(i-2)/2 


4+(i-l) 


4+(i-l)/2 


Constraints 


3+(i-l) 


3+(i-l) 


3+(i-l) 


3+(i-l) 


Solution 


Unique 


NO 


Undetermined 


NO 



TABLE IV: The number of parameters on which LTB models depend expanding the luminosity dis- 
tance to i-th order are shown in the smooth case without age constraint and according to the presence 
of a homogeneous bang function tb{r). 



4- Analysis of the relevant cases 

The information about the existence of a solution of the inversion problem can be obtained 
by comparing the number of the free parameters to the number of constraints, but some cases 
require some further analysis. 

5. C* models 

The distinction between odd and even order of the Taylor expansion for the case comes 
from the fact that if i is odd, then we need to expand to {i — l)-th order in r, which is an even 
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Smooth with age constraint 


Matching to i-th order in Dl{z) 


Even 


Odd 


LTB type 


Smooth, tb ^0 


Smooth, ti, = 


Smooth, tb ^0 


Smooth, tb = 


Free parameters 


4+(i-2) 


4+(i-2)/2 


4+(i-l) 


4+(i-l)/2 


Constraints 


3+(i-l)+l 


3+(i-l)+l 


3+(i-l)+l 


3+(i-l)+l 


Solution 


NO 


NO 


NO 


NO 



TABLE V: The number of parameters on which LTB models depend expanding the luminosity distance 
to i-th order are shown in the smooth case with age constraint and according to the presence of a 
homogeneous bang function tb{r). 



number, can be included in the expansion, and so we have some new free parameter. On the 
contrary if i is even we should expand to the (i — l)-th in r, which is odd, and because of the 
smoothness conditions the (i — l)th order cannot be included, but we still have an additional 
constraint coming from the D^^^z) matching condition, making the inversion problem more 
difficult to solve. Intuitively this can be interpreted as the fact that smooth models have less 
degrees of freedom to satisfy the necessary observational constraints, and going to higher order 
in the redshift expansion does not always add new free parameters because of the smoothness 
conditions which require all the odd derivatives for the the functions k{r) and tb{r) to vanish at 
the center. These are the general properties of the solution of the inversion problem for Dl{z) 
for models, and the differences respect to the case: 

• If we include the age constraint there is no solution. This is due to the same reason why 
there is no solution also in the case studied in the previous section, i.e. D2 cannot be 
matched. 

• If we do not include the age constraint and tb{r) ^ there is a slightly undetermined 
solution when i is odd , but only one parameter is left free, the one which is not fixed by 
the age constraint. This residual freedom is less than in model, where it increases as 
with the matching order i. In the case the matching order is even the solution is unique. 

• If the age constraint is not included and tb{r) — there is no solution, while in the 
case there was a unique solution. 
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IX. CONCLUSION 



We have extended previous studies about central smoothness in LTB models, by considering 
the implications of higher order smoothness conditions. Using a low redshift expansion we have 
analyzed in detail what are the the constraints that a LTB model has to satisfy in order to 
be in agreement with Hq, luminosity distance and the age of the Universe observations. The 
constraint on the age of Universe is actually exact, since the analytical solution can be used to 
impose it. We have also given a systematic analysis of the existence of solutions to the inversion 
problem in the different cases depending on the central smoothness and the presence of a bang 
function t;,(r), comparing the solutions investigated in previous studies to the more general 

case. For both and models there is not solution because of the impossibility to match 
the second order D2term in the expansion of Dl{z). One difference is in the case in which the 
age constrain is not included and the bang function is zero, in which there is a unique solution 
for models but no solution for the C* case. Another difference is in the case in which the age 
constrain is not included with non zero bang function, in which the solution is undetermined for 
both and models, but the latter ones have much less residual parametric freedom. Our 
formulation of the inversion problem uses dimensionless parameters to define the LTB model, 
which make evident the freedom of fitting any Hq observed value, correcting some previous 
claims that this last observable would be enough to rule out these models. The dimensionless 
formulation of the inversion problem we have derived has the advantage to express the LTB 
solution directly in terms of observable quantities, which is more transparent than deducing 
Hq from other parameters of the void model as done in [s^ for example. These results can be 
used to define a general algorithm to look for LTB models in agreement with both low and high 
redshift observations, using H^^'^ as the fundamental scale in terms of which to define all the 
other parameters and functions determining the model. Our results imply that any LTB model 
able to fit luminosity distance data, satisfy the age constraint and fit some other observable is 
either not mimicking exactly the AC DM red-shift space observables theoretical predictions or 
it is not C°° smooth. 

While the present paper was focused on pressureless LTB models, the same arguments could 
be applied to LTB models in presence of a cosmological constant A as proposed in 43l-l46| . and 
in this case even small inhomogeneities could have an important effect on the apparent value A 
and on the equation of state of dark energy. We will consider this important effects in separate 
works, underlying nevertheless their general importance since they do not require large voids 
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and should be taken properly into account even in standard ACDM models. 
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